Low temperature spin fluctuations in geometrically frustrated Yb3Ga5 0i2. 
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In the garnet structure compound Yb3Ga50i2, the Yb^^ ions (ground state effective spin S' — 
1/2) are situated on two interpenetrating corner sharing triangular sublattices such that frustrated 
magnetic interactions are possible. Previous specific heat measurements have evidenced the devel- 
opment of short range magnetic correlations below ~ 0.5K and a A-transition at 0.054K (Filippi et 
al. J. Phys. C: Solid State Physics 13 (1980) 1277). From ^™Yb Mossbauer spectroscopy measure- 
ments down to 36 mK, we find there is no static magnetic order at temperatures below that of the 
A-transition. Below ^ 0.3 K, the fluctuation frequency of the short range correlated Yb"^"*" moments 
progressively slows down and as T ^ 0, it tends to a quasi-saturated value of 3 x 10^ s~^. We 
also examined the Yb''+ paramagnetic relaxation rates up to 300 K using Yb perturbed angular 
correlation measurements: they evidence phonon driven processes. 



PACS numbers: 76.80.-|-y, 75.50.Lk 
I. INTRODUCTION 

For most crystallographically ordered compounds con- 
taining magnetic ions, the limiting low temperature mag- 
netic ground state involves long range magnetic order 
where the spin fluctuations die out as T — > 0. For some 
particular lattice structures however, the geometric ar- 
rangement of the magnetic ions is such that it may not 
be possible to simultaneously minimise all pairs of inter- 
action energies. The resulting frustration may then lead 
to a situation where long range order does not occur [1,2] 
and where the presence of a large number of low energy 
states leads to the continued presence of spin fluctuations 
as T ^ 0. 

Systems with frustrated interactions that are of current 
interest include the kagome lattice [3,4], where the ions 
are arranged on a motif of corner sharing triangles, the 
pyrochlore lattice [5] , where the ions are arranged on cor- 
ner sharing tetrahedra and the garnet lattice (R3T5O12) 
[6], where the rare earths (R) form two interpenetrat- 
ing, non-coplanar, corner sharing triangular sublattices. 
This geometry does allow frustration to be operative pro- 
vided there is a suitable combination of the nature and 
the size of the rare earth anisotropy and of the sign of 
the interionic interactions. A number of the rare earth 
garnets appear to evidence a conventional long range or- 
dered state [7] suggesting that in these cases, frustration 
plays a negligible role. In fact, to date, frustration has 
been reported to play a major role in only one garnet, 
Gd3Ga50i2 [8-13], where the S-state Gd^+ ion has a 
very small intrinsic anisotropy and where the dominant 
coupling is antiferromagnetic. 

Amongst the garnets made with the non S-state 
rare earths, Yb3Ga50i2 is unusual in that the Yb'^+ 



ground state shows only a relatively modest crystal field 
anisotropy (see below). Specific heat data has evidenced 
a broad peak centred near 0.2 K, attributed to short 
range correlations and a sharp A-peak at 0.054 K, initially 
attributed to the onset of long range magnetic order [14]. 
We have carried out ^^°Yb Mossbauer spectroscopy mea- 
surements down to 0.036 K in order to examine the be- 
haviour of the Yb'^"*' moments as the temperature is low- 
ered through that of the A-transition and to examine the 
low temperature spin dynamics. We also report ^''^Yb 
perturbed angular correlation measurements, carried out 
from 14 to 300 K, which provide information concerning 
the thermal dependence of the Yb'^^ fluctuation rates in 
the paramagnetic region. 



II. BACKGROUND PROPERTIES 

The single phase polycrystalline sample was prepared 
by heating the constituent oxides to 1100°C four times 
with intermediate grindings. 

In the garnet lattice (space group Ia3d), the rare earth 
site point symmetry is orthorhombic [mmm). The crys- 
tal field acts on the Yb'^+, ^^7/2 state to leave a ground 
state Kramers doublet which is very well isolated from 
the excited Kramers doublets [15]. For Yb"^"*" ions di- 
luted in Y3Ga50i2, the ground doublet g-values are gx 
= 3.73, gy = 3.60 and gz = 2.85 [16,17], and the wave 
function is derived from the cubic Fy state (g = 3.43). In 
Yb3Ga50i2, the Yb'^''" g-values should be quite similar. 
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FIG. 1. Specific heat in Yb3Ga50i2 reproduced from Ref. 
[14]. 

A magnetic 4f-shell contribution to the specific heat 
is visible below ^ 0.5 K (see Fig.l): there is a broad 
peak centered near 0.2 K, followed by a A-anomaly at 
0.054 K. The total electronic entropy change below IK 
is very close to the value Rln2 expected for an isolated 
Kramers doublet, and only about 10% of it is released 
at the A-transition [14]. If we attribute the broad peak 
to exchange driven short range correlations, then the ex- 
change energy scale is ~ 0.2 K. The susceptibility data 
[14,18] follow a Curie- Weiss behaviour down to ~ 1.0 K 
with a small paramagnetic Curie- Weiss temperature [Op 
= 0.05 K) corresponding to a net ferromagnetic inter- 
action. Below ~ 1.0 K, in the region where the broad 
specific heat peak occurs and where the magnetic corre- 
lations develop, the thermal dependence of the suscepti- 
bility falls below that corresponding to the extrapolated 
Curie- Weiss dependence. This behaviour evidences the 
presence of magnetic correlations which are antiferomag- 
netic. In Yb3Ga50i2, there is thus evidence for the exis- 
tence of both ferromagnetic and antiferromagnetic inter- 
actions. The small value for Op suggests the two types of 
interactions have comparable strengths. The magnetic 
frustration in Yb3Ga50i2 that is evidenced in this re- 
port, thus appears to be linked to the presence of antifer- 
romagnetic correlations within a Heisenberg-like system 
on triangular sublattices and to the presence of inter- 
actions with opposite signs. The isomorphous compound 
Gd3Ga50i2 where frustration is also operative, evidences 
a dominant nearest neighbour interaction which is an- 
tiferromagnetic and other interactions with competing 
signs [8]. 



III. "°YB MOSSBAUER MEASUREMENTS 

A. General features 

The ™Yh Mossbauer absorption measurements (Ig = 
0, le = 2, = 84keV, Icm/s corresponds to 680 MHz) 
were made down to 0.036 K in a ^He-'*He dilution refrig- 
erator using a neutron activated TmBi2 source displaced 
with a triangular velocity sweep. Selected spectra at 4.2, 
0.15, 0.075 and 0.036 K are shown in Fig.2. The last 
two temperatures are situated either side of that of the 
specific heat A-transition (0.054 K). 
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FIG. 2. i™Yb''*+ Mossbauer absorption in Yb3Ga50i2. At 
4.2 K, the fitted line was obtained using a relaxation lineshape 
appropriate for paramagnetic fluctuations. At 0.15, 0.075 and 
0.036 K, the data fits were obtained using a lineshape appro- 
priate for hyperfine field fluctuations. As the temperature 
decreases, the absorption broadens and its centre of gravity 
moves towards velocities which are more negative than the 
isomer shift value shown by the dashed line. 

At 4.2 K, the absorption takes the form of a broad 
Lorentzian shaped line, with a half-width at half max- 
imum, r = 4mm/s, significantly broader than the exper- 
imental halfwidth of the TmBi2 source, Fq = 1.35 mm/s. 
As the temperature is lowered to ~ 0.25 K, the absorp- 
tion line remains Lorentzian shaped and keeps essentially 
the same width, but it moves slightly towards negative 
Doppler velocities. This means the centre of gravity 
of the absorption no longer corresponds to the value of 
the isomer shift, which is close to 0.0 mm/s relative to 
the TmBi2 source. Below ~ 0.3 K, the absorption line 
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moves markedly towards more negative velocities, it also 
progressively broadens and becomes slightly asymmetric 
(spectra at 0.075 and 0.036 K in Fig.2). No resolved hy- 
perfine structure is visible at any temperature, even at 
temperatures below that of the A-transition. This indi- 
cates there is no "static" long or short range magnetic 
order and that magnetic fluctuations persist down to the 
lowest temperatures. In relation to the characteristic fre- 
quency scale of the present ^™Yb Mossbauer measure- 
ments, the absence of a well resolved hypcrfinc structure 
means that the fluctuation frequency of the Yb^+ mag- 
netic moments remains above the threshhold value of ~ 
3xl0®s~^. The quantitative analysis of the fluctuation 
rate is presented in the next section. 

B. Quantitative analysis 

The broad, single-line-nature of the absorption in the 
paramagnetic region at 4.2 K arises because the Yb'^"'" 
magnetic hyperfine splitting is "motional narrowed" by 
the fast fluctuations of the Yb'^"'" magnetic moment. It 
is possible to extract the fluctuation frequency from the 
measured lincshape by using a paramagnetic spin relax- 
ation model based on a perturbative approach [19], pro- 
vided the magnetic hypcrfinc tensor A is known. With 
the approximation of local axial symmetry, wc obtained 
the components of this tensor from a ^^*^Yb Mossbauer 
measurement on Yb'^^ ions diluted into Y3Ga50i2, 
where because the dilution removes the spin-spin cou- 
pling, the fully resolved hyperfine splitting is observ- 
able. We obtained the values: Az/h = 738 MHz and 
A±/h = 952 MHz, corresponding respectively to g- values: 
g2=2.82 and g^=3.63. These values are essentially equiv- 
alent to those previously measured by electron spin reso- 
nance (gz=2.85 and g±={gx+gy)/'2=3M [18]). We note 
that with an isotropic A tensor and in the fast relaxation 
limit (i.e. when A/h <C l/r, where l/r is the Yb^"*" 
paramagnetic spin relaxation rate), this lineshape model 
leads to a single line having a Lorentzian shape with a 
dynamical half broadening given by: 

ATu = l{A/nrr. (1) 

On fitting the data at 4.2 K using the perturbation 
approach [19] and the axially symmetric A tensor com- 
ponents, we obtain the Yb''+ paramagnetic spin fluctu- 
ation rate: l/r ~ 3.8 x 10^° s~^. This rate is constant 
between ~ 0.25 K and the highest measurement temper- 
ature of 80 K, indicating that the driving mechanism is 
the temperature independent spin-spin coupling between 
the Yb^+ ions. A rough estimate of the strength of 
this coupling, using the relation: U/t ^ ksT^x, yields: 
Tex ~ 0.3 K, consistent with the temperature of the 
broad maximum of the exchange induced specific heat 
(~ 0.2K). 



The shift of the centre of gravity of the absorption away 
from the temperature independent isomer shift value ev- 
idenced in Fig.2 is analogous to that previously observed 
at low temperatures in YbAlOa [20] and in YbBcia [21]. 
It seems to be related to the inadequacy, in these cases, 
of the perturbative relaxation lineshape when the driving 
mechanism is the exchange spin-spin interaction. As dis- 
cussed in Ref. [21], it is likely that the matrix elements 
of the relaxation operator have non-vanishing imaginary 
parts [22,23], and these generate lineshapes whose spec- 
tral signature is an anomalous shift of the centre of grav- 
ity of the absorption. The fact that the anomalous shift 
increases as the temperature decreases could be linked to 
the growing influence of the spin-spin correlations, hence 
to the growing inadequacy of the standard perturbative 
relaxation lineshape. The relaxation rates obtained from 
the perturbation analysis are not influenced by the shift 
of the centre of gravity of the absorption since the relax- 
ation rate is linked to the real part of the matrix elements, 
i.e. to the width of the absorption line. 

At very low temperatures (T < O.IK), the shape of 
the experimental absorption can no longer be correctly 
reproduced by the paramagnetic relaxation model. Since 
the specific heat data (Fig.l) show that magnetic corre- 
lations are present in this temperature region, wc fitted 
the experimental data below 1 K using a relaxation model 
involving hyperfine field fluctuations which were consid- 
ered in the random phase approximation (RPA) [24]. A 
hyperfine field {Hhf) is indeed present at the -^^"Yb nu- 
cleus when the Yb^+ moments are short (or long) range 
correlated and the size of the field is proportional to that 
of the Yb'^^ moment. The time dependent Hamiltonian 
for this relaxation lineshape is: 

H^Hq- gnl^uHhf J2 (2) 

j=l,N 

where Hq is the quadrupolar hyperfine Hamiltonian, £f„ 

the gyromagnetic factor of the excited nuclear state, /i„ 
the nuclear Bohr magneton, the summation is over the 
N directions among which the hyperfine field fluctuates 
and fj (t) is a random function of time with appropriate 
values corresponding to the different possible forms of 
the Hamiltonian. For Hq , we take the (very small) value 
we have determined from -'^^'^Yb measurements on Yb'^+ 
substituted into Y3Ga50i2. Then, the lineshape depends 
on a single dynamic parameter, the fluctuation frequency 
of the hypcrfinc field (l/r)/iy, and on the choice of the 
directions between which the fiuctuations take place. We 
obtain very satisfactory fits to the data below 0.2 K (solid 
lines in Fig.2 at 0.036, 0.075 and 0.15K) by assuming the 
hyperfine field fluctuates between the three principal di- 
rections of the local coordinate frame, that is, between 
OX, OY and OZ, the principal directions of the electric 
field gradient tensor. In the Hamiltonian (2), this corre- 
sponds to A''=6 and j = ±X, ±Y and ± Z. If we assume 
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the hypcrfinc field fluctuates along only one of these axes, 
we obtain much poorer data fits. In the rapid relaxation 
rate limit, this model yields a line of Lorentzian shape 
with a dynamical half-width: 



' aMkf 



(3) 



The 0.036 K lineshape is well broadened and this en- 
ables us to obtain both the magnitude of the hyperfine 
field and its fluctuation rate. We flnd: H^f c:^ 140(10) T, 
which corresponds to a Yb^+ moment of ~ 1.4 (for 
i70Yb3+, l^B yields a hyperfine field of 102 T) and 
l/r)/j/ = 3 X 10^ s~^ The value for the Yb^+ moment is 
not far from the mean value expected both from the aver- 
age g-tensor (3.4) and from the saturated magnetisation 
measured at 0.09K in Ref. [14] (1.7/Ub). In the fits of 
the spectra up to 0.2 K, we used an intrinsic half-width: 
ro=1.35mm/s and we assumed the fluctuating hyperfine 
field has a size which remains constant at the value 140 T 
derived at 0.036 K. This assumption should hold as long 
as the correlations are well developed, but it cannot be 
ascertained that it is completely correct up to 0.2 K. 
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FIG. 3. Thermal variations, in Yb3GaBOi2, of the Yb^+ 
hyperfine field fluctuation frequency extracted from the 
^^"Vb Mossbauer spectra. The dashed line is the law: 
a(l/r)ft/ = O.SksT. 

The thermal variation of {1/T)hf is shown on Fig. 3. As 
the temperature is lowered over the range 0.2 to O.IK, 
the frequency decreases approximately linearly with a 
law: h{l/T)hf = O.^ksT, and then below about O.IK, 
it tends to saturate towards the value 3 x 10^ s~^. There 
is essentially no difference between the rates either side 
of the specific heat A-transition (0.054 K). The T-lincar 
dependence of (1/t)/j/ only pertains to a limited temper- 
ature range, and it must be kept in mind that it somehow 
depends on the validity of the assumption about the con- 
stant magnitude of the fiuctuating hyperfine field. We 
note that a linear variation has also been encountered 



theoretically for the case of a frustrated Hciscnbcrg py- 
rochlore antiferromagnet [5]. The decrease in the relax- 
ation rate followed by a saturation is similar to the be- 
haviour observed (by fiSK [13]) in the isomorphous com- 
pound Gd3Ga50i2. 



IV. i^^YB PERTURBED ANGULAR 
CORRELATION (PAC) MEASUREMENTS 

The PAC measurements provide the Yb^+ param- 
agnetic relaxation rates in the temperature range 14 
- 300 K. The measurements were made using the 91- 
1094 keV 7 — 7 cascade from the ^"^^Lu ^''^Yb /3 decay. 
The i^^Lu nuclei are obtained by proton irradiating the 
sample, which is then annealed at 800° C to remove ir- 
radiation defects. The intermediate level of the cascade, 
namely the 1172 keV rmclear excited level of ^"^^Yb with 
spin / = 3 and half-life 8.3 ns, is used to observe the per- 
turbation of the 7 — 7 directional correlations due to the 
hyperfine interactions. In the case of a static quadrupo- 
lar or magnetic hyperfine interaction, oscillations are ob- 
served in the time evolution of the perturbation factor 
R{t) [25]. 
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FIG. 4. "^Yb PAC spectra at 14 K and 300 K, fitted to an 
exponential decay. The first channel, spoilt due to prompt 
coincidences, was removed from the fits. 

Up to ^300 K, we observe the perturbation factor R{t) 
does not show oscillations but it is rather an exponential 
function of time: R{t) = Aexp{—nt) (see Fig. 4). This 
suggests a dynamic hyperfine interaction. Paramagnetic 
relaxation within the ground Yb'^+ doublet, with effective 
spin 1/2, leads to such an exponential decay, in the fast 
relaxation limit: <C l/r, where ^^^^ is the hy- 

perfine constant of the intermediate / = 3 level of ^^^Yb 
and l/r the electronic spin fluctuation frequency. Then, 
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in a manner analogous to Eqn.l, the damping rate fj, is 
given by [26]: 

M=^(Ai^V?^)V. (4) 

The value for the hyperfine constant of the 7 = 3 
level of ^''^Yb can be obtained by scaling the mean 
A™ value appropriate for the Yb'^"'" ground state 
{A^™/h ~ 884 MHz) with the nuclear g-factors. We ob- 
tain: A^'^'^/h ^ 565 MHz. The thermal variation of 1/r 
derived from that of n is shown in Fig. 5 which also shows 
the ^'^°Yb Mossbauer derived values at 4.2 and 80 K. The 
two sets of values agree quite well. The relaxation rate 
is constant up to ~ 150 K, then increases monotonically. 
The observed thermal dependence can be fitted to the 
sum of a temperature independent spin-spin term and 
an exponential term associated with a two-phonon real 
process through the excited crystal field states (Orbach 
process) [27]: 

^ = (iu + i^exp(-^), (5) 

where A is the energy of an excited crystal field level. 
The fit yields: (I/t)^^ = 3.8 x lO^'^ s'^ B = 3.3 x lO^^ 
s~^ and A = 880(50) K. This last value is in very good 
agreement with the mean energy distance between the 
Yb3+ ground doublet and the three closely spaced excited 
crystal field doublets (~ 850 K) [15]. Additional PAC 
measurements concerning the thermal variation of the 
electric field gradient at the ^^^Yb nucleus above 300 K 
are reported in Ref. [28]. 
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FIG. 5. Thermal variation of the Yb 4f shell magnetic 
fluctuation rate in Yb3Ga50i2 derived from perturbed angu- 
lar correlations measurements (black dots). Two values ob- 
tained from the ^'^"Yb Mossbauer spectra at 4.2 and 80 K arc 
also shown (open squares). The solid line is a fit to a sum of 
spin-spin and spin-phonon driven relaxation rates (see text). 



V. DISCUSSION AND CONCLUSIONS 

Previous specific heat measurements [14] have shown 
that, in Yb3Ga5 0i2, magnetic correlations develop be- 
low ~ 0.5 K and that a A-transition occurs at 0.054 K. 
Our present investigations, using ^™Yb Mossbauer spec- 
troscopy down to 36 inK and ^^^Yb perturbed angular 
correlations measurements up to 300 K, provide insight 
into the correlations and the dynamic behaviour of the 
Yb"^"*" spins over a wide temperature range. 

Above ~ 0.5 K, the Yb^+ magnetic moments undergo 
paramagnetic fluctuations. Up to ~ 150 K, the fluc- 
tuation rate has a temperature independent value of 
~ 3.8 X 10^° s~^ and the driving mechanism is the ex- 
change interaction between the the Yb^"*" spins. Above 
150 K, additional temperature dependent relaxation oc- 
curs through coupling to phonons according to a two- 
phonon Orbach process involving the excited crystal field 
states near 850 K. 

The low temperature Yb'^"'" magnetic correlations show 
up in the ^^°Yb Mossbauer measurements below ~ 0.3 K 
through changes in the lineshape. The fluctuation fre- 
quency of the correlated moments decreases as the tem- 
perature is lowered and below 0.1 K, it tends to a quasi- 
saturated value of 3xl0^s~^. This is a quite high value 
and in fact, Yb3Ga50i2 is the only known compound 
where the T — > spin fluctuation rate is rapid enough to 
fall within the ^^"Yb Mossbauer spectroscopy frequency 
window. The T ^ state in Yb3Ga50i2 is therefore 
a dynamic short range correlated spin-liquid state. On 
crossing the temperatures of the specific heat A-anomaly, 
there is no significant change in the ^™Yh Mossbauer 
lineshape. Usually, a transition to a long range mag- 
netically ordered phase reveals itself in the Mossbauer 
spectra by the appearance of a well defined magnetic hy- 
perfine splitting. This does not appear in the present 
case, clearly showing there is no long range order at tem- 
peratures below that of the specific heat peak. This peak 
has the intriguing characteristic that the associated en- 
tropy gain is very small, i.e. it amoimts to about 10% 
of the total Rln2 entropy gain associated with the Yb"^+ 
ground state Kramers doublet. 

Some examples are already known of frustrated sys- 
tems where the specific heat peak has an associated low 
entropy. In the pyrochlore compound Gd2Ti207, there is 
an entropy gain of 50% of i? In 8 at the transition at 1 K 
[29], which has been shown by neutron diffraction to in- 
volve long range magnetic order [30]. Yb2Ti207 presents 
an entropy gain of about 20% of i?ln2 [31] at the tran- 
sition at 0.25 K which has been shown to be associated 
with a first order change in the fluctuation rate of the cor- 
related Yb''+ moments [32]. In Yb3Ga50i2, the reduced 
size specific heat peak is linked neither to the appearance 
of long range order nor to a mesurable change in the fluc- 
tuation rate of the short range correlated moments. In 



5 



fact, we observe a slowing down of the fluctuations of the 
correlated moments below ~ 0.3 K, i.e. at temperatures 
well above that of the specific heat peak, and the fluctu- 
ation rate is quasi temperature independent in the region 
where the specific heat peak occurs. 

Usually for frustrated systems presenting a phase tran- 
sition, the temperature at which this transition occurs 
is sizeably lower than the temperature associated with 
the strength of the interionic interaction. The energy 
scale of this interaction may be estimated in different 
ways: from the paramagnetic Curie- Weiss temperature, 
from the temperature where the specific heat evidences a 
magnetic contribution (in the present case, the tempera- 
ture of the broad maximum) and from the paramagnetic 
spin-spin relaxation rate. In Yb3Ga50i2, the last two 
methods lead to an interaction of equivalent strength 0.2 - 
0.3 K, whereas the paramagnetic Curie temperature (Op), 
is smaller (0.05 K). The correspondence between 9p and 
the temperature of the specific heat peak appears to be a 
mere coincidence. The precise origin of this low entropy 
specific heat peak within a spin liquid phase remains an 
unresolved issTic. 
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